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Abstract 

We describe the possible scenarios for the evolution of a thin spherically symmet¬ 
ric self-gravitating phantom shell around the Schwarzschild black hole. The general 
equations describing the motion of the shell with a general form of equation of state 
are derived and analyzed. The different types of space-time R±- and T^-regions 
and shell motion are classified depending on the parameters of the problem. It is 
shown that in the case of a positive shell mass there exist three scenarios for the 
shell evolution with an infinite motion and two distinctive types of collapse. Analo¬ 
gous scenarios were classified for the case of a negative shell mass. In particular this 
classification shows that it is impossible for the physical observer to detect the fan- 
tom energy flow. We shortly discuss the importance of our results for astrophysical 
applications. 
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1 Introduction 


Recent observations of both the type la high-redshift Supernovas (SNs) and 
Cosmic Microwave Background (CMB) anisotropy at the small angular scales 
strongly indicate in favor of the acceleration of the universe expansion at the 
present epoch [1], The simplest possibility for the acceleration expansion of 
the universe is an existence of the cosmological constant A with an equation of 
state w = p/p = — 1 [2], The cosmological constant provides the satisfactory 
explanation of the cosmic dynamics but encounters the fine tuning problem. 
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An alternative explanation is the existence of dark energy in the form of 
a specific scalar held (quintessence) whose equation of state is varying with 
time (see e. g. [3,4]). Such a model allows to construct the so called ‘tracker’ or 
‘attractor’ cosmological solutions which resolve in particular the cosmological 
fine tuning problem [4]. 

One of the peculiar feature of the dark energy is a possibility of the phantom 
energy equation of state p + p < 0. This phantom energy violates the weak 
energy dominance condition. In the case of phantom energy the cosmological 
scenario with the ‘Big Rip’ is possible when cosmological phantom energy 
density grows at large times and disrupts finally all bounded objects up to 
subnuclear scale [5]. The other peculiarity is the diminishing of the black hole 
mass due to phantom energy accretion [6]. The phantom energy is usually 
associated with the phantom or ghost fields: e. g. the scalar fields with a wrong 
sign kinetic term [5,7]. The phantom energy has some peculiar properties 
in the framework of QFT in the curved space-time [8]. The thermodynamic 
properties of the phantom energy is also rather unusual [9]. 

The existence of phantom energy is not excluded by the nowadays observa¬ 
tions. The measurements [10] of the distances and host extinctions of the 
230 SN la provide the constraints on the dark energy equation of state, 
— 1.48 < w < —0.72. In [11] the data set containing 172 type la supernovas 
are analyzed in the model independent manner and it was shown that the 
presence of the phantom energy with —1.2 < w < —1 is preferable for the 
recent epoch. The Chandra telescope observations [12] of the hot gas in the 26 
X-ray luminous dynamically relaxed galaxy clusters provides w = — 1.20to;28, 
which is also in favor of the phantom energy. 

The evolution of dark energy are considered usually in relation to the cosmo¬ 
logical problems. However the local evolution of self-gravitating dark energy 
may be quite different from the cosmological one. This is because of the non¬ 
linearity of the General Relativity equations. The three dimensional analytical 
treatment is possible only in very restrictive cases, e. g. in the case of the sta¬ 
tionary accretion onto black hole of the dark energy considered as test fluid, 
i. e. with a negligible self-gravitation [6]. self-gravitating fluid one must go 
to some simplified models. One of the analytically treatable approach with 
a fluid self-gravitation is taking into account is the thin shells model. The 
theory of thin shells in General Relativity was developed by W. Israel [13] 
and developed then by many authors (see e. g. [14] for review and references). 
The problem of thin shell analysis is greatly simplified in the case of spherical 
symmetry. The aim of this paper is to consider several scenarios for the thin 
spherically symmetric phantom shell evolution. 

The paper is organized as follows. In Sec. 2 the general concepts of spherically 
symmetric gravity are outlined with special attention to Schwarzschild space- 
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time. In Sec. 3 the specific equations of motion for thin shells are obtained. 
In Sec. 4 the evolution of shell with phantom equation of state is analyzed 
and different types of motion are classified. In Sec. 5 we briefly discuss the 
obtained results. Throughout the paper we use the units h — c — 1. 


2 General Theory 

2.1 Spherically symmetric gravity 


A spherically symmetric manifold is a direct product of a two-dimensional 
manifold M 2 and two-dimensional sphere S 2 , that is, M 4 = M 2 x 62 • The line 
element of any spherically symmetric space-time can always be written in the 
form 


ds 2 = g a pdx a dx 13 = Adt 2 + 2 Hdtdq + Bdq 2 + R 2 (t, q)dSl 2 (1) 


with the signature (+, —, —, —). Here t and q are correspondingly the timelike 
and spacelikc coordinates, A, H and B are functions of t and q only, and 
R{t , q) is the radius of a two dimensional sphere (in the sense that the area of 
the sphere equals to 4ttR 2 ), 

dti 2 = d6 2 + sin 0 2 d(f) 2 ( 2 ) 


being the line element of the unit sphere. For the given space-time the co¬ 
efficients A, H and B are not uniquely defined. One can transform the line 
element ( 1 ) to the new coordinate system which conserves explicitly the spher¬ 
ically symmetric form of the metric: 

t = t(t,q), q = q{t,q). (3) 


Unlike the metric coefficients in M 2 , the radius R(t, q ,) is invariant under the 
transformation (3). The other very important invariant is constructed from 
the partial derivatives of R as follows 

A = 7 a ' 3 R,aRj3, (4) 

where is inverse to the two-dimensional metric tensor 7 a p. This invariant 
is nothing more but the square of the normal vector to the surface R = const. 

If we know two invariant functions R(t, q) and A (t, q) we know the line element 
of the spherically symmetric space-time up to the gauge transformations and, 
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therefore, its local structure. To construct the global manifold we need some 
additional principle. Physics provides us with it. From the point of view of 
physicists any space-time should be geodesically complete [2], that is, every 
timelike and null geodesics should start and end either at infinities or at the 
singularities. 

The function A (t, q) brings a nontrivial information about a space-time struc¬ 
ture. Indeed, in the flat Minkowskian space-time A = —1, all the surfaces 
R = const are timelike and therefore, R can be chosen as spatial coordinate 
q — R on the whole manifold. But in the curved space-time A is no more 
constant and can in general be both positive and negative. The region with 
A < 0 is called the /2-region, and the radius can be chosen as a radial coor¬ 
dinate q. In the region with A > 0 the surfaces R = const are spacelike (the 
normal vector is timelike), and the radius R can be chosen as a time coor¬ 
dinate t. Such regions are called the T-regions. The R- and T-regions were 
introduced by Igor Novikov. But this is not the whole story. It is easy to show 
that we can not get R = 0 (“dot” means a time derivative) in a T-region. 
Hence it must be either R > 0 (such a region of inevitable expansion is called 
T + -region) or R < 0 (inevitable contraction, a T_-region). The same holds 
for R- regions. They are divided in two classes, those with R' > 0 (“prime” 
stands for a spatial derivative) which are called i? + -regions, and TL-regions 
with R' < 0. These, R- and T-regions are separated by the surfaces A = 0 
which are called the apparent horizons. The apparent horizons can be null, 
timelike or spacelike. 

Thus, the curved spherically symmetric space-times may in general have a 
rather complex structure, a set of R±- and T±-regions separated by apparent 
horizons A = 0. In the next subsection we consider one of the important 
example of spherically symmetric manifolds. 


2.2 Schwarzschild space-time 


The solutions to the vacuum Einstein equations consist of only one-parameter 
family. In the curvature coordinates (q = R) the metric of a Schwarzschild 
space-time has the form 

ds 2 = F(R)dt 2 - F-\R)dR 2 - i?W, (5) 

where 

F(R) = 1 - m > 0 (6) 

it 
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Fig. 1. Carter-Penrose diagram for Schwarzschild space-time. Every point represents 
a sphere. Here /_ and J+ mean past and future null infinities respectively. The past 
and future timelike and spacelike infinities are i_ and iq respectively. Regions R 
and T are separated by the two apparent horizons (the future and past ones) which 
are null surfaces. In our case the future horizon coincides with the event horizon 
defined as the first null geodesics which does not reach infinity. Its time reversal is 
called a particle horizon. The T_-region is called the black hole. The T + -region is 
called the white hole, and the A_-region is called the wormhole respectively. 

and G is the Newton gravitational constant. This static metric describes 
a space-time outside a spherically symmetric body with a mass m called 
also a Schwarzschild mass. The metric (5) has a coordinate singularity at 
R = 2 Gm. It is related to the static nature of the line element (5) and im¬ 
possibility to synchronize clocks of the observers who are static at the spatial 
infinity (R, = const —> oo) and those (who are nonstatic) in the region 
R < 2 Gm. Moreover, it appears that the manifold described by the line in¬ 
terval (5) is not geodesically complete. It can be shown that the maximally 
extended Schwarzschild manifold consists of four parts. Furthermore, it is 
possible to choose such coordinates, which put infinities at finite distances. 
Schwarzschild geometry in such coordinates is represented by the Carter- 
Penrose diagram shown in the Fig. 1. An another useful representation of 
the Schwarzschild space-time is the so called embedding diagrams. One can 
consider Schwarzschild metric with t = const and 0 = 7r/2. It is easy to show 
that this is the metric on a hyperboloid embedded in the three dimensional 
flat space 

r = sk + 2Gm - (7) 

So, an embedding looks like it is shown in the Fig. 2 In the following we 
will use schematic version of embedding which is called embedding diagram 
(right panel in the Fig. 2). As one can see from the Fig. 2 the wormhole is 
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Fig. 2. Schwarzschild spacetime at fixed time and 8 and its embedding diagram. 
The throat width is proportional to Schwarzschild mass. 

separated from the R + -region (or the black hole exterior) by the throat called 
also the Einstein-Rosen bridge. Of course, the above definitions are by no 
means general, but they are sufficient for our purposes. (The general definitions 
require powerful mathematical tools and enormous number of predefinitions.) 
With the inclusion of matter sources the Schwarzschild solution is valid only 
outside their boundaries. In the case of the complete Schwarzschild manifold 
the sources are considered as concentrated at the past and future singularities 
at R = 0. The manifold is called the eternal black hole. 


3 Thin shells 


One of the most important feature of General Relativity is that the equa¬ 
tions of motion of matter fields are incorporated into the Einstein equations. 
The Einstein equations of General Relativity are nonlinear partial differen¬ 
tial equations. This means that the motion of test particles or fields on the 
given background will in general be different from that of the matter for the 
self-consistent solutions. It makes analysis very complicated. To obtain some 
definite results we have to choose the simplest possible model. This is a self- 
gravitating thin shell. In this section we derive equations of motion for thin 
shell. 

Let us now introduce the notion of thin shell. Because we will be dealing 
only with the timelike spherically symmetric thin shells we adjust the very 
nice generally covariant formalism derived by W. Israel [13] to the case of 
interest (see, e.g. [14]). Let us choose some hypersurface E divided the whole 
space time into two parts, “in” and “out”. With this hypersurface E can be 
connected some special coordinate system called Gauss normal coordinates. 
In our simple case the line element written in these coordinates takes the form 

ds 2 = dr 2 — dn 2 — f? 2 (r, n)dRt 2 . (8) 
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where r is the proper time of the observer sitting fixed on E, the coordinate 
n grows from the “in” to the “out”-region in the outer normal direction to 
the hypersurface E, R(r,n) is the radius of the sphere (in the sense that a 
sphere area equals to AnR 2 ) and dRl 2 is the line element of the unit sphere, 
dU 2 = d6 2 + sin 2 Odxfi 2 . The hypersurface is situated at n = 0, and 

= dr 2 — p 2 (r)dh2 2 , ^ 

p( T ) = R(t, 0) 


The hypersurface E is called the singular shell if some energy momentum 
tensor is concentrated on it, namely Tf = SfS(n) + ..., Sf is the surface 
energy-momentum tensor of the shell (i, k= 0,2,3), otherwise the hypersurface 
is nonsingular, in our case due to the spherical symmetry the only nonzero 
components of Sf are Sq and S 2 = S^. 

The introduced earlier invariant A equals 

A = R 2 -R, 2 n (10) 


and 


R 


in 


E — 



in) 


where a = +1 if radii increase in the direction of the outer normal, and a — — 1 
if radii decrease. Evidently, a — +1 in Z? + -region and a = — 1 in i?„-region. 
Thus on the equation of motion of the shell a can change sign only in T-region, 
or in the region where its motion is forbidden. 


The subsequent procedure is very simple. Keeping in mind that the metric 
itself is continuous but some of its derivatives could undergo a jump across 
the shell, we integrate Einstein equations and obtain (the nontrivial result is 
only for the (q) and (|) components) after some algebra 

—\P 2 - A ta - - A out = 8 irGpSS, (12) 

P P 


2 ( 71 , 


Vp^a- 




A 


out 


a/p 2 ~ a i 


-P 


, CTout p+ 
\J P 2 —Aout 

(1 + A out ) + 


_ ( 1 _i_ A. ^ _ _ gout^ 

2 P\J p 2 -A in m 2p^/ p 2 - A out 

+4vrG'p( (out )T^ T") = 8 tt GS$. 


(13) 
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The continuity equation for the energy-momentum tensor is transformed to 

^ + ^(Sq 0 - S%) + < out) T 0 n - (in) T n _ q ( 14 ) 

The third equation is a differential consequence of the first two ones. But very 
often it is convenient to use all three of them. 

In what follows we will consider the thin shell in vacuum. So, both inside and 
outside the shell we will have the Schwarzschild metric with different masses 
(because the shell has its own mass which is added to the inner mass) and our 
equations become: 


^"in \j -^in ^"out \j “I - -^out 4:7tG pS^ ^ 


(15) 

” a / C 0\2 , q 2/^i2 n0 q 2 G( m m + m out) 

p— 4 7T Ct p(6 0 ) + 87T Ct pb Q b 2 

Am 2 S'! 

(16) 

87t 2 p 5 (Sq) 3 

ss + 2 As§ - Sj) = 0, 


(17) 


where A = — F = — 1 + AAy We wrote the second equation in a somewhat 
different (twice squared) form, which is more suitable for us. The information 
about a global geometry (signs of cr in and cr out ) is already contained in the 
equation (15). 


4 Evolution of phantom shell 

Now let us apply the above theory to phantom shells. Consider a simple (but 
a rather general!) linear equation-of-state relation S q = kS 2 - In the phantom 
case k > 1. The solution of (17) is 

S° = Cp 2(fc_1) , (18) 

where we call the constant C the ’shell power’. Denote also 

x = (19) 

After some manipulations with (15) and taking into account (19) one gets the 
following two equations: 

. 2 1 G (. Sm 2 

P = -H- {m in + m OVLt ) + — -b x 

p \ 4x 
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( 20 ) 




Fig. 3. Schematic curves p 2 (p) (up) and p(p) (down) according to (20) and (21). 
The case C 2 > Cq and C 2 < Cq corresponds to the left and right panel respectively. 
The region p 2 < 0 is forbidden for motion. 


P=~: 


(min + m ou t) H-2(2 k - l)x 

x 


( 21 ) 


and the sign conditions 

din = sign 5 m + 8tt 2 Gp 3 (Sq) , (22) 

o- out = sign 5m - 87r 2 Gp 3 (S° ) (23) 

or in a more convenient form 

<r in = sign (5m + 2x ), (24) 

(T out = sign (5m - 2x). (25) 


Here, we denote 5m = m out — m in . 

The curves p 2 (p) and p(p) are shown in the Fig. 3. There are possible several 
evolution scenarios for the shell. When the curve p 2 (p) lies completely above 
the p axis, there is only the infinite motion. On the other case, when the 
function p 2 (p) has roots, it is possible both the finite motion and infinite one. 

As one can see from equation (19), it is convenient to define the parameter 
space of the problem using m in , m out , k > 1/4 and p as free parameters. Then, 


9 



instead of search for conditions imposed on p, it is enough to find conditions for 
the parameter C. This conditions will define the global geometry. Now, let us 
construct these conditions. First of all, the change of sign of the acceleration p 
in (21) occurs when p — 0. This corresponds to the quadratic equation, whose 
positive root is 


x 0 


m in + m out + yf ( 4 k - 1 ) 2 5m 2 + 4 m in m out 
4(2 k - 1) 


(26) 


The corresponding value of p according to (19), denoted by po- Consider the 
sign of p 2 (p 0 ). Let us introduce the parameter C 0 that p 2 (p 0 ) > 0 when 

C 2 > C 2 (27) 


and p 2 (po) < 0 when 


C 2 < C 2 


(28) 


The explicit value for Cq is 


1-4fc 

} ’ ■ (29) 


x\(4k-l)G 


(4 k - 1) (m in + m out ) + yf (4 k -1) 2 5m 2 + 4m in m out 


Consider at first the case of 8m > 0 . According to (24) there must be always 
<T in = +1- The <T out changes sign at x = X\ = 5m/ 2 and <r out = — 1 if x > X\. 
Thus, in the case x —> oo, one has cr ou t = — 1, and, if x — > 0 , then cr out = +1. 
The value of p corresponding to x\ is denoted by p\. From (21) it is easy to 
see that 

■pix ,) = < 0 . ( 30 ) 

Pi 


From this follows the important conclusion that p\ < po- This can also be 
proved by the direct comparison of Xq and x i. From (20) one obtains 


P 2 {x i) 


^ ^ 2 Gm out 
Pi 


(31) 


10 




c- 


c, 


Fig. 4. Fraction (Ci/Co ) 2 as a function of p and k. 

Let us denote by C\ such a value of C that p\ = 2Gm out at C — C\\ 

8m 


Cf = 


2 2(2fc+l) 7r 2 m 4fc-1 ' 


(32) 


In the case pi > 2Grn out one has 

C 2 < C 2 (33) 

and cr out changes sign in the forbidden for motion part of f?-region. And vice 
versa, if pi < 2Gm out , then 

C 2 > C 2 (34) 


and cr out changes sign in T-region. Now we have to define which of the con¬ 
ditions Cq > C\ or C{ < Gq is true. This can be done by considering limits 
for the parameters in (29) and(32) or by direct numerical calculation of the 
fraction (Ci/Co) 2 as a function of p — Srn/rri out and k. The results of such 
calculation is shown in the Fig. 4. It is clear now that Cf < Cq. 

It is obvious from the above analysis, that there exist three possible scenarios 
for the shell evolution in the case 8m > 0. In short, these scenarios are the 
following: 
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out 



Fig. 5. Infinite inflation or collapse (symmetric in time). The left throat is the more 
narrow then the right one because 5m > 0. 

i’+ r=0 r=0 *+ i’+ r=0 i+=i+ r=0 i + 




Fig. 6. The Carter-Penrose diagrams for a collapsing (left) and inflating (right) shell. 




Fig. 7. The Carter-Penrose diagrams for a collapsing shell (left) and for a shell 
moving from infinity to infinity (right). 

(1) Infinite motion. The sign of cr ou t changes in the T region. The shell power 
obeys the following inequality: 


C‘ 


>cl 


(35) 


The embedding diagram for this scenario is shown in the Fig. 5. The 
Carter-Penrose diagrams for the case of collapse and inflation are the 
same up to the time reverse (see Fig. 6). 

(2) There exist the turn points and cr out changes sign in T-region: 

Cl < C 2 < Cl (36) 

In this case the embedding diagrams are the same as in the previous case. 
The corresponding Carter-Penrose diagrams are shown in the Fig. 7 

(3) There exist the turn points and <7 out changes the sign in a part of P-region, 
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Fig. 8. The left diagram represents the 
left turn point. In this case the shell is 
shell goes from the past infinity to the 
right from the right turn point. 



case when the shell is at the left from the 
collapsed finally. At the right diagram the 
future infinity when it evolves toward the 


i’+ r=0 i’ + =i + r=0 i + 



Fig. 9. The Carter-Penrose diagrams for the collapsing shell (left) and for the shell 
moving from infinity to infinity (right). In the both cases the shell evolves in T±- 
and i?_-regions. 


which is forbidden for the motion. 


C 2 < Cl (37) 

In this case we have two embedding diagrams (see the Fig. 8). The Carter- 
Penrose diagrams are shown in the Fig. 9 As follows from the above 
analysis if 5m > 0, the right diagram in the Fig. 9 is the only case when 
the shell shows itself in the i? + -region. 

Consider now the 5m < 0 case. The ’density’ Sq is assumed to be always 
positive. So the negativity of 5m is caused exclusively by the gravitational 
mass defect. The cr in changes sign at x = x 2 = —5m/(2), and <j in = +1 if 
x > x 2 . At the same time cr out = — 1 Thus, we can conclude that in the case 
x —> oo, one has cr in = +1 and cr in = ^1 if x —> 0. Denote according to (19) 
the corresponding to x 2 value of p by p 2 . From (21) it is easy to see that 

p(x 2 ) = < o. (38) 

P 2 
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in 


out 


Fig. 10. Infinite inflation or collapse (symmetric in time). The right throat is the 
more narrow then the left one because 5m > 0. 

From (20) one obtains 


P 2 (x 2 ) 


^ 2Gm m 
P2 


(39) 


Let us denote also by C 2 such a value of C that p 2 = 2 Gm in at C = C 2 : 




5 m 

2 2(2»+l) x 2 m 4J t -l- 


(40) 


In the case p 2 > 2Grn m one has 


C 2 < Cf 


(41) 


and the considered region is i?-region. And vice verse, if p 2 < 2Gm m then 


C 2 > cl 


(42) 


and one has T-region here. Just as in the case 5m > 0, the fraction (C 2 /Co) 2 < 
1 because the 2d-graph is analog of the Fig. 4 with the replacement /i —> —fi. 
As a rsult, we can conclude that in the case 5m < 0 there exist two evolution 
scenarios: 

(1) Infinite motion. The sign of cri n changes in T region. The inequality for 
C'-parameters are 

C 2 > Cl (43) 

The embedding diagram for this scenario is shown in the Fig. 10. The 
Carter-Penrose diagrams for the case of the collapse and inflation are the 
same up to the time reverse, (see the Fig. 11). 
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Fig. 11. The left Carter-Penrose diagram represents a collapsing shell. The right one 
represents an inflating shell. 



Fig. 12. The left Carter-Penrose diagram represents collapse of the shell. At the 
right one the shell moves from infinity to infinity. In both cases shell evolves in T±- 
and R --regions (for a distant observer). 

(2) There exist the turn points and cr out changes sign in T-region. 

C 2 < C 2 < C 2 . (44) 

In this case the embedding diagrams are the same as in the previous case. 
The corresponding Carter-Penrose diagrams are shown in the Fig. 12 

(3) There exist the turn points and cr out changes sign in a part of i?-region, 
which is forbidden for the motion. 

C 2 < C 2 . (45) 

For this case we have two embedding diagrams (see the Fig. 13). The 
Carter-Penrose diagrams are shown in the Fig. 14. 

We see again that for the all scenarios in the case of 5m < 0 the shell evolves 
under horizons and cannot reach a distant observer living in the i? + -region. 


5 Conclusion 


Discussion 

We considered a dynamics of phantom thin shell surrounded Schwarzschild 
black hole. The motivation for this work is the fact that in many physically 
interesting situations in cosmology and astrophysics the essential role was 
played the full account for gravitational backreaction. In our case of phantom 
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in 


out 




out 


in 


Fig. 13. The left diagram represents situation when the shell is on the left from the 
left turn point. In this case the shell is collapsed. At the right diagram shell goes 
from past infinity to future infinity when it evolves to the right from the right turn 
point. 


i’ + r=0 r=0 i + i’+ r ~0 r=0 i + 



Fig. 14. The left Carter-Penrose diagram represents collapse of the shell. At the 
right one shell moves from infinity to infinity. In both cases shell evolves in T±- and 
i?_-regions 

shells such a backreaction may appear crucial for formation of the global ge¬ 
ometry of the space-time. The matter is that in General Relativity any type 
of energy is gravitating. That is, not only energy density but also the tension 
and pressure are gravitating. The pressure plays a twofold role. The positive 
pressure causes both repulsion and attraction, the later is due to its contribu¬ 
tion to the gravitating source. The negative pressure, on the contrary, leads to 
the gravitational repulsion (the famous example is the deSitter space-time). 
Hence the phantom shell is even more repulsive. And indeed, we show that 
the global geometry of the system consisting of the Schwarzschild surrounded 
by the phantom shell is the wormhole-like type in all but one cases. In the 
wormhole-like type geometry the distant observer cannot see the shell at all, 
they are separated by the throat (Einstein-Rosen bridge). The only exception 
is the case of the bound motion with 5m > 0. But, though the distant observer 
may see the shell it can not register the energy flux of the shell. 

We are sure that despite the very simple character of our model the result 
obtained should be taken into account in doing calculation in cosmology and 
astrophysics when phantom energy is present. 
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